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ABSTRACT 

We investigate aspects of the cosmological evolution of FR II radio galax- 
ies, focusing first on the abilities of models to match data for linear-sizes, radio 
powers, redshifts and spectral indices. Here we consider modifications to the 
theoretical models we had treated earlier, primarily by accounting for the growth 
of the radius of hotspots with source size. Better fits to the distributions of most 
of the data in three low-frequency surveys can be found with sensible choices 
of model parameters but no model yet considered gives a good match to all of 
the survey data simultaneously, nor does any do a good job of producing the 
spectral index distributions. The observational datasets are too small to com- 
pletely discriminate among the models. We calculate the volume fraction of the 
"relevant universe" cumulatively occupied by the expanding radio galaxy lobes 
over the quasar era, when these powerful radio galaxies were much more com- 
mon, and when they have been argued to play an important role in triggering 
galaxy formation and spreading magnetic fields and metals. We found the cumu- 
lative relevant volume filling factor of radio galaxies to be ~ 5%, so we conclude 
that these impacts are smaller than previously estimated but that they are still 
significant. 

Subject headings: galaxies: active — galaxies: luminosity function, mass function 
— large-scale structure of universe — methods: statistical — radio continuum: 
galaxies — surveys 
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Introduction 



Classical double radio galaxies (RGs) with extended lobes, or the Fanaroff-Riley Class II 
(FR II) sources, constitute the more powerful RG population. The quasar era (i.e., be- 
tween redshifts ~ 1.5 and 3) is distinguishable by rises in the comovine; luminosity densities 



of RGs. Quasars and other Dowerful AGN (DunloD & 


^eacockl 199d J 


ackson k Wall 1999: 


Willott et al. 2001: Ueda et al. 2003: 


Grimes et al. 2004 


Hopkins e 


t al. 


2006h, as well as high 


star and galaxy formation rates (e.g., 


Lillv et al. 19961: Madau et al. 


199J 


V. Bouwens & Illineworth 


2006: Sawicki &; Thompson 2006a b 


1. These similar evolutions suggest that powerful RGs 



can have substantial impacts on the for mation, distribution and evo 



large scale structures of the universe (e .g., Gopal-Krishna fc Wiital 



2001 



Rawlings fc Jarvis 



2004 



Gopal-Krishna fc Wiital l2003al. GK W03a: iGopal-Krishna fc Wiital l2003bl. GKW 03b 



2001 



ution of g alaxies and 



GKWOL lKronberg et al. 



Gopal-Krishna. Wiita. fc Ostermanll2003l. G KWO: lGopal-Krishna. Wiita. fc Barail booi GKWB04; 



Levine fc GnedinI l2005l : ISilkl 12005 ). Observational inyestigations 



show r en iarkable radio-optical alignment in high- 2: RGs (e.g 



1988a 



bl:lDev et al. lll997l:IVen emans et al.l[20o3 . 12OO5I : ICreve et al 



McCarthv et al.lll987 : Chambers et al. 



2006 



Overzier et al. 



Zheng et alibood : kiiki et alibood l Some studies attribute this to FR II RGs triggering 



2006 



ex- 



tensive star forraation in a rnulti-phase intergalac tic medium (IGM) (e.g.. iBegelman fc CiofB 



19891 : iDe Youndll989l : lReeslll989l : iChokshil 119971 : GKWOl). In addition, the expanding RG 



lobes could easily have infused significant magnetic fields into the IGM (G KWOl, GKWO, 



GKWB). Si milar conclusions w' e re dra wn from different lines of argument by lKronberg et al. 
( I2OOII ) and iFurlanetto fc Loebl (1200 ll ). These radio sources born through the quasar era 



could also have contributed toward metal enhancement of their environments (GKW03b 



GKWB04), since observations fe.g.. [Dietrich et al. 



20041 : IShaplev et al.ll2004l : IPieri et al. 



2006 



2003 



Prochaska et al 



2006 



Schave et al.l l2003l: lAguirre et al. 



Tripp et al.ll2006l ) require 



an efficient mechanism for spreading metals widely at early cosmic epochs. 

A comprehensive study of the impact of RGs on various events in the cosmological 
history of the universe requires reliable quantitative estimates of the relevant volume filling 
fraction of RGs and their active lifetimes. A prerequisite for a more accurate computation of 
this volume impacted by radio lobes is a good model of the evolutio n of radio sources, for both 
individual sources and as a function of redshift. In a recent work iRawlings fc Jarvis! (120041 ) 
agreed that RG lobes will penetrate much of the relevant universe, but they argued that this 
may often shut off star formation by expelling gas fr om prqtoclu sters. However, they assume 



a single phase medium, unlike our picture or that of iRees 



is not surprising. Many recent obseryational studies (e.g.. 



Jl989h. so this negative conclusion 



Bohringer et al. 1995: Birzan et al. 



20041 : IChoi et al.ll2004l : lO'Dea et al.ll2004l : iRevnolds et al.ll2005r ) show depressions in X-ray 
surface brightness coinciding with radio lobes, cavities and buoyant bubbles in clusters of 
galaxies; these are clear signatures of the interactions of radio sources with their surrounding 



- 3 - 



hot intracluster gas on scales of several tens to hundreds of kiloparsecs from the AGN. 



Our aim is to probe in more detail the cosmological impact of RGs on the growth of 
structures. To do so, we must first develop an improved but essentially analytical model 
for the evolution of FR II RGs as they age individually and as their number s vary with 
cosmological epoch. As our first step toward this goal we (IBarai fc Wiitall2006l . BW) com- 
pared three fairly sophisticated analy tical models for the dynamical and lobe pow er evo- 
lution of FR II RGs, those g i ven b y iKaiser. Dennett-Thorpe, fc Alexandeii (119971 . KDA) , 
Blundell. Rawhnes. fc WillottI Jl999l . BRW) and lManolakoufc Kirkl J2OO2I . MK): the reader 
should have some familiarity with these papers or see iBarai for more details. The 

source linear-size evolution in the BRW and MK models essentially follow the KDA prescrip- 
tion. However the models differ in the ways the relativistic particles are injected from the 
jet to the lobe, and in treatments of loss terms and particle transport. So there are some 
significant differences in their predictions for observed powers (P) as functions of source size 
(D) and redshift (z). 

The [-P^-D] evolutiona ry tracks of model radio sources have been used in KDA, MK, and 



m 



Machalski et al.l (j2004al Jbl) as the main way to look for consistency between observational 



data a nd RG evolution models. The radio sky simulation prescription in lBlundell. Rawlings. &: Willott 
( 119991 . BRW) adds new dimensions to the observed parame ter space. Deriving a RG redshift 
distribution function from the work of lWillott et al.l (I2OOII ) on the radio luminosity function 
(RLF), BRW prescribed a virtual RG survey technique. This involves generation of a huge 
initial population of sources over extended cosmic epoch according to pre-defined distribution 
functions in redshift, jet power, source age, and orientation. These large numbers of sources 
then evolve through their individual lives where they interact with their environments and 
undergo energy losses. Only a few simulated sources are detected in the virtual surveys 
when mathematical flux limits, corresponding to observations, are ir nposed. We adopted ob- 
serva tional samples from the redshift-comp lete subsam ples: 3CRR (ILaing. Rilev. fc Longair 
1983h . 6CE jRawhngs. Bales, fc LacvlboOlh and 7CRS (IMcGilchrist et al.lll99oh . of the Cam- 
bridge fiux-limited radio catalogs; see Table 2 of BW. The models, relevant observations and 
the multi-dimensional Monte Carlo simulation methodology leading t o the virtual surveys 
are discussed in detail in §§2 - 4 of BW, and given in more detail in iBarail (120061 ). Using 
the virtual surveys and any radio lobe power evolution model (KDA, BRW or MK), one can 
get P, D, z, and spectral index, a (with the convention Pi, oc z/~°), values for the simulated 
model radio sources detected in the pseudo-surveys. The distributions of the characteristics 
of these simulated RGs can then be compared to observational data to test the model. In 
order to assess the success of a theoretical model, we perform several statistical tests and 
present the quantified results for the various model fits. In BW, our main conclusion is that 
none of the heretofore published models provides an adequate description of the data. 
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As our next major step toward the goal of isolating good RG models, we have varied 
the radio lobe power evolution models in the present work. We consider modifications to the 
original models of BRW and MK, by incorporating a variable hotspot size growing with the 
source age. Here, we perform analogous simulation-based virtual surveys on the modified 
models, and present the corresponding statistical results . In O we describe the models: 
Modified BRW (MBRW), Modified MK (MMK), K2000 JKaiseilT200ol ) and the alternative 
radio luminosity function (RLF). We give the statistical test results of the model simulations 
in ^ which we discuss in ^ 

One important goal of this work has been to address the question of what fraction of 
the relevant volume of the universe (the volume containing most of the cosmic baryons) 
did the radio lobes occupy during the quasar era (GKWOl). The warm/hot intergalactic 
medium (WHIM) comprises the main repositor y of cosmic baryons w hich can potentially 
collapse to form star clusters or galaxies (e.g., ICen fc Ostrikerl 119991 ). So the radio lobes 
need to penetrate a significant portion of this "relevant volume of the universe" occupied 
by the WHIM filaments in order to have a significant role in impacting star formation and 
spreading magnetic fields and metals (e.g., GKWB04; BGKOW04 and references therein). In 
^we calculate the total volume filled by the RGs over the quasar era (when their population 
peaked) as a fraction of the relevant volume of the universe. Our conclusions are in ^ 



2. Modified Radio Galaxy Evolution Models 

The major modification to the BRW and MK models involved allowing the hotspots to 
grow in size as a source ages and expands . The data used to make a sensible modification 



are taken from I Jeyakumar fc Saikial (120001 . hereafter JSOO), who studied the dependence of 
sizes of hotspots on overall source sizes for a sample of FR II sources which included both 
compact steep spectrum and larger-sized sources spanning a projected source size range from 
about 50 pc to nearly 1 Mpc (Fig. [1]). We parameterize the hotspot radius, 

rhs = rhso + f{L), (1) 

where Vhso is some normalizing initial hotspot radius, and f{L) is a power law expression of 
the total linear size L of the source. We chose rhsg such that the hotspot of a source grows to 
= 2.5 kpc when the total linear size is L = 200 kpc, since these are reasonable averages 
of the actual values and = 2.5 kpc was the constant value assumed by BRW. 

The hotspot and source angular size data are adopted from JSOO. We follow these 
authors and find the average angular hotspot size for each source. This is the geometric 
mean of each hotspot (major and minor axes) sizes and, for those sources with hotspots 
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Fig. 1. — Quadratic cu rve fit to the hotspot size vs. linear size, [Dhs^L] data of 
Jevakumar fc Saikial (§2]). 



detected on botli sides, tlie aritlimetic average of their sizes. We convert the angular hotspot 
sizes and separation s to corresponding projected linear sizes using the consensus cosmology 
( jSpergel et al.ll2006l ). Total linear sizes (L) are obtained by assuming an average angle to 
the line of sight of 39.5° (following KDA, as done for the [P-D] tracks in BW §5.1). The 
hotspots are much smaller than the total source sizes and are assumed to be spherical, so 
projection effects are negligible for them. 

We performed various least-squares fits to the \og{rhs) vs. log(L) data: a single straight 
line; a single quadratic; two straight lines with a break at 20 kpc; two straight lines with a 
break at 1 kpc. Although all these fits are satisfactory, a quadratic fit to the data gave the 
lowest reduced and is shown in Fig. [H So the hotspot size was taken to be growing with 
the source size as 

\og{Dhs) = Chs, + yi log L + y2 (log Lf . (2) 

Then the hotspot radius was given by Eq. ([T]) using f{L) = Dhs/2. The best fit values 
of the coeffic ients are Ch,,, = —3.199, yi = 1.053, and ?/2 = —0.0306, which gave a reduced 
Xr = 0.2720. iPerucho fc Marti] (120031 ) presented a dynamical model for compact-symmetric 
and FR II sources where the hotspot size gr ows self-similarly with linear size, and can get a 
fit which is in agreement with the results of iJeyakumar fc Saikial (120001 ). 



In the modified BRW (MBRW) model the hotspot size (and area, A/, 



vrr 



rises 



according to Eqs. ([T]) and ([2]). It otherwise follows the prescription of the original BRW 
model summarized in §3.3 of BW. During a source's evolution several additional quantities 
(which were fixed for a source when a constant hotspot size was assumed) also varied with 
its age. These include the hotspot pressure, phs, the hotspot magnetic field, Bhs, the "slow" 
and "fast" break frequencies, and Ubf, and the corresponding critical Lorentz factors, jbs 



- 6 - 



and 7b/. 



The modified MK (MMK) model also incorporates the same rising hotspot size (and 
area). In the MK model, the characteristic time, to, (Eq. 6 of MK) when the size of the 
head was comparable to the hotspot size, depends on "an initial" hotspot area Ahs{to), which 
we must now distinguish from the normal rising hotspot area Ahs- We found this "initial" 
hotspot area using r/is(to) = 0.02 kpc. We chose this value as it gave the best results when 
we compared the statistics of 6 MMK simulation runs done using rhs{to) = 0.01 — 0.06 kpc, 
computed at intervals of 0.01 kpc. So in the MMK model to is 



3cl-'^cA 



1 1/« 



/is (to) 



(r, + 1) (5 - 



Qo 



3/(4+/3) 



(3) 



with A 



hs{to) 



hs{to) 



TT (0.02 kpc)^, po is the central density of the ambient gas, Oq is its 
scale length, [3 is its radial density index, Qq is the jet power, is the adiabatic index of 
the external environment, a is (4 + /?) / (5 — /3), and ci is a model dependent (but weakly 
varying) constant (see BW for details). The MMK model otherwise follows the prescription 
of the original MK model as described in §3.4 of BW. 



Kaiserl (120001 . K2000) proposed a modification to the KDA (IKaiser. Dennett-Thorpe, fc Alexander 
19971 ) model, whi ch essentially follows the fo rmulation in KDA but which has the following 
differences. The iKaiser fc Alexander! (119971 ) and the KDA model considered a cylindrical 
geometry for the cocoon, where the hotspot pressure ph drives the source expansion along 
the jet axis and the expansion perpendicular to the axis is governed by the cocoon pressure 
Pc- Hence it was assumed in the KDA model that 



Pc 



(4) 



where Rt is the axial ratio, or the ratio of the length of the source to the full width of a lobe 
half way down t he jet. In subsequent studie s of the shocked gas flow between the bow shock 
and the cocoon. IKaiser &: Alexanderl (119991 ) empirically fitted {ph/Pc) as functions of (3 and 
Rt- Such investigations appeared to show that the ratio in Eq. (jl]) was an overestimate. To 
correct this problem K2000 claimed that a better empirical formula was given by 



^ = (2.14 

Pc 



0.52(3) R'^''-'-''^. 



(5) 



The source birth function over redshift. 



p{z) oc exp 



1 / z- Zq 

2 V cj,. 



(6) 
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is the common gaussian functional form of all the Radio Luminosity Functions (RLFs) we 
considered. A peak redshift of Zq = 2.2, and a standard deviation of = 0.6, following 



Willott et al.l (120011 ). was adopted in the model simulations we performed in BW and in 



most of this paper. To explore the effect o f different assumed redshift d is tribut ions we also 
considered extensively the RLE given by iGrimes. Rawlings. fc WillottI (120041 ). which has 
zo = 1.68 4 and = 0.447, us ing the two-population generalised luminosity function from 
Table 5 of ICrimes et~aD J2004h . 



3. Model Results 
3.1. [P-D] Tracks 

The power (P) vs. linear-size (D), or, [P-D] tracks of the MBRW, MMK, and K2000 
models, as well as those of the original BRW, MK, and KDA models are shown in Fig. [2l 
The modified models follow the same general trends as do the original models described in 
§5.1 of BW. The tracks are generated using the modified models (described in ^ with the 
default values of parameters for dynamical and power evolution from each of the original 
models (given in Table 1 of BW). Each source was evolved at frequency u = 151 MHz. For 
this figure the total linear sizes were converted to the projected sizes assuming an average 
viewing angle to the line of sight of 39.5° (following KDA). 

The rates of steepening of the tracks are significantly different in the three new models. 
The MBRW track is less steep than the original BRW track. Among the three original 
models, KDA, BRW and MK, BRW gave the worst fit to the data, when compared with 
respect to K-S statistical tests (BW). If the reason for this can be identified with the fact 
that BRW gave the steepest [P-D] tracks (Fig. 1 of BW), then this "shallowing" of the 
tracks in the MBRW model would imply that the K-S statistical fits to the data should be 
better for the MBRW model (which is indeed true). The MMK track is slightly steeper than 
the original MK track, implying less of a difference between their fits to the data. The K2000 
track is much flatter than the corresponding typical tracks of any other model. 



The "youth-redshift degeneracy" (e.g.. iBlundell fc Rawlingslll999l ) is clear in these [P- 



D] tracks. A high-power, high- redshift, source shows a faster fall off in its specific 151 MHz 
luminosity with time, and can even fall below the limiting flux of a radio survey at a younger 
age, as compared to a lower-power source at a lower z. 
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Fig. 2. — [P-D] tracks of a fiducial source having jet power, Qq = 1.3 x 10^^ Watts, and 
at redshift, z = 0.5. The curves of different hnestyles correspond to the tracks predicted by 
different models as labeled. The crosses on each track denote source lifetimes of 1, 10, 20, 
30, 90, 100, 150, 200, 250, 300 Myr. 

3.2. Preliminary Statistics: 1-D K-S Tests 

We used 1-dimensional Kolmogorov-Smirnov (1-D K-S) statistics as a first quantitative 
test. Each of the distributions of key characteristics [P,D,z,a] of the radio sources detected 
in the simulated surveys were compared to those of the sources in the real radio surveys 3C, 
6C, and 7C. The K-S probabilities, V, that the two data sets being compared are drawn 
from the same distribution function, were taken to be figures of merit of each simulation. To 
quantify the overall success of a model, we added the K-S probability statistic for comparisons 
of P, D, z, a (i.e., V{P) +V{D) -\-V{z) +V{a)) for the three surveys, weighting the statistic 
of a survey by the square-root of the number of simulated sources detected in that survey, 
and denote this overall figure of merit as V[p^D,z,a]- The second figure of merit we employ, 
denoted P[p,2Z),2,a], adds the K-S statistic probabilities for P and z to twice those for D, i.e. 
V{P) + 2V{D) + V{z) + V{a) using the same weighting method. More details can be found 
in §5.2.1 of BW. For each model some parameters were chosen as "better" (in providing 1-D 
K-S fits to the data); for these, further simulations and additional statistical tests ( §3.31 and 
13. 4p were done. 

The procedures followed in doing the simulations and in presenting the results are briefly 
discussed below; for complete details see §5 of BW. Table [T] gives our results for the MBRW 
model; Table [2] gives those for the MMK model, and Table [3] those for the K2000 model. 
The tables follow the same format and pattern as the corresponding tables 3, 4 and 5 of BW 
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for the original models, and should be compared to them. All the statistical test (K-S and 
correlations) results for the modified models include a 1 kpc cutoff of source size, i.e., the 
statistics are calculated by excluding sources with linear size, D < 1 kpc (see §5.2.3 of BW). 

An initial ensemble, generated using the default parameters from BRW for the RG 
population generation, was evolved according to each of the modified models (MBRW, MMK 
and K2000). The sources in the simulated surveys (produced according to the prescription 
of §4.3 of BW) were compared to the data samples of the 3C, 6C and 7C catalogs. We 
examined the 1-D K-S test statistics of the first entry (the very first 3 rows) of Tables [H [2] 
and [3] (modified model results) and compared those to the first entries of Tables 4, 5 and 
3 of BW, respectively (original model results). From this single comparison we find that 
the MBRW model is significantly better (the combined 1-D K-S probabilities are ~ 10 — 20 
times higher, which is ~ 4a better) than the original BRW model, while the default MMK 
model is only slightly better than the original MK model. The K2000 model produces much 
worse K-S fits than did the original KDA model, and for the reasons discussed in §l]we did 
not explore the K2000 variation any further. 

In search of further improvements of the combined 1-D K-S statistics, we varied the beam 
power distribution function of the sources generated in the initial population by checking 
steeper exponents, x, in the power-law of the initial jet power distribution, where, p{Qo) oc 
Qq^ between Qmin and Qmax- For the MBRW model the overall statistics improved the most 
for X = 3.0, while for the MMK model x = 2.6 gave better fits. These values of x were then 
used for the later simulations. 

The initial population generated with x = 3 (but otherwise using the BRW prescription), 
was evolved according to the MBRW power evolution model. The corresponding 1-D K-S 
statistics are given in the third entry in Table [T], and show the improved fit (compared to 
the first entry of that table). To search for possible further improvements we varied the 
other parameters prescribing the power evolution in the models. Simulated surveys were 
constructed using the parameter listing given in Table [1] (each variation done one at a time) 
of the MBRW power evolution model. The total 1-D K-S statistics, V[p^D,z,a] and V[p^2D,z,a], 
as seen from Table [T], are comparable to or better than the original BRW model results 
(Table 4 of BW). 

In order to find the best-fit maximum RG age, we performed simulation runs using 
initial populations with x = 3 (for MBRW), and x = 2.6 (for MMK), and then varied 
TMax- The Tjifax that gave the highest mean values of the combined statistics (V[p^D,z,a] and 
V[p^2D,z,a]) was chosen as the best maximum age at the considered x. In the MBRW model 
the highest mean statistics were seen at Tmgx = 300 Myr. The MMK model performed its 
best at TMax = 150 Myr. Hence we used initial populations with the above "optimal" x and 
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Tuax pairs for each model in subsequent runs. The 1-D K-S results of a subset of the first 
set of runs employing these "optimal" initial populations and ages but varying some of the 
power evolution model parameters to alternate values are in Tables [T] and [2] for the MBRW 
and MMK models, respectively. 

Upon examining these preliminary results, only those cases that gave any improvement 
in statistics over the default case or were essentially as good as the default were considered 
further. For these parameter sets three more runs (making a total of four) were done using the 
same large population but with different pseudo-random seeds, yielding different ensembles 
of simluated samples. Then the means and standard deviations of the relevant 1-D K-S 
statistics were found. Some "2-change" cases, i.e., models where two "superior" parameter 
variations (those giving high 1-D K-S probabilities) were simultaneously employed also were 
explored. As seen from the tables showing the individual 1-D K-S statistic probabilities 
(Tables [Hand [2]), often several of the 12 K-S probabilities for some cases give acceptable fits, 
and a few would give very good fits, but it is difficult to find a single model where all are 
good fits. In other words, the modified models also do not give good simultaneous fits to the 
[P, a] data from all three of the radio surveys considered (3C, 6C and 7C). 

The 1-D K-S st atistical test results for the simulation runs of different models using the 



Grimes et al.l (120041 ) RLF are given in Table HI 



3.3. 2-D K-S Tests 

We performed additional statistical analyses on the modified models, so as to make 
more robust comparisons both between them and with the original models. This was done 
in a fashion similar to that of §5.3 of BW for the original BRW, KDA and MK models. 
For each modified model the 1-D K-S best-fit parameter variation cases, i.e., those which 
gave the highest combined probability 'P[p,D,2,a]; according to the results from the previous 
sub-section were selected. 

The 2-D K-S test results for both the default versions of the modified models and the 
parameter sets (denoted as "varied") giving the highest total 1-D K-S probability for each 
model, are given in Table [51 The results are listed in a similar way as are the 1-D K-S 
statistics in previous tables, and are discussed below ( §4.ip . 
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3.4. Correlation Coefficients 

Spearman partial rank correlation coefficients were calculated for those cases for which 
the 2-D K-S tests were done. We combined the [P, D, z, a] data from the 3 surveys: 3C, 6C 
and 7C III, for the actual observations and for the model simulations, and computed correla- 
tions between them (see §5.3.2 of BW for the limitation to 7C 111). This was mainly done in 
order to subdue the tight [P-z] correlation present in a single fiux-limited complete survey, 
and to thereby discover any correlations which exist between the other source characteristics. 

TableOgives the four-variable Spearman partial rank correlation coefficients {rp£,,za-, ^Pz,d, 
etc.) which were computed on the combined data and on the simulated results from the 
modified models. We also examined the corresponding 2-variable and 3-variable correla- 
tions. Probably the most important result is that the 2-variable correlation, rp^, is always 
negative; however, when the 4- variable correlation between P and D is found, (i.e., with 
the effects of z and a removed), a small positive correlation is seen between P and D (i.e., 
^PD,za > 0). These are the same as trends seen in the [-P--D] correlations for the original 
models. 



4. Discussion of Simulation Results 

The K2000 model produces very flat [P-D] tracks (Fig. [1]) and the 1-D K-S statistical 
results obtained in the multi-dimensional Monte Carlo simulations (Table[3]), are all very poor 
compared to any of the other models studied here or in BW. So we conclude that the K2000 
model cannot well reproduce the trends of observed properties in the low frequency radio 
surveys. The reason for such poor behavior of this model appears to be that the K2000 model 
was designed specifically to describe the environments and ages of three local, and rather 
atypical, FR 11 sources. As this model is biased toward describing special environments, the 
parameters used in K2000 cannot really be applied globally and hence it is not surprising 
that this model cannot explain the cosmological evolution of RGs. Since this aspect is the 
crux of our investigations, we do not consider the K2000 model any further. 

The primary modification made to the other models, the incorporation of a growing 
hotspot size, produces the following major results. The MBRW model (results in §3.2[ 
Tabled]) is a substantially better ffi to the data than the original BRW model (§5.2 in BW), 
as the total 1-D K-S probability is better by ~ 4cr in the default case, and by ~ 2.5o" in 
the "best-ffi" case, of MBRW when compared to BRW. The 1-D K-S proba bilities for a 



are sometimes better (especially for 7C) and in a few cases (see appendix of iBarail l2006l ) 



approach the value 0.01 for which a model spectral index ffi is not ffimly rejected. The 
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MMK model produced fitting statistics ( §3.2[ Table [2]) which are better than or comparable 
to the original MK model fits (§5.2 in BW). 

We have explored the modified models through our extensive multi-dimensional Monte 
Carlo simulation procedures and parameter variations in the modified models. But similarly 
to what we found for the original models in §5 of BW, we found that no modified model 
gives acceptable fits to all the source characteristics, [P, D, z, a], for all the three surveys 3C, 
6C and 7C, simultaneously. 

Steepening the power law index for the initial beam power distribution (Eq. 2 of BW) 
to a; = 3 (from x = 2.6 used by BRW) while using the default maximum age Tmgx = 500 
Myr, improved the 1-D K-S statistics for the MBRW model, as can be seen from Table [H 
However, in the MMK model x = 2.6 gave better results, as compared to x = 3 (Table [2]). 
This is the only model we examined which gave better 1-D K-S fits with the beam power 
distribution index set to a; = 2.6. Simulations done by co-varying Tmqx and x gave better 
fits at a maximum age of 300 Myr (MBRW) and 150 Myr (MMK), when combined with the 
above "optimal" values of x. These "best-fit" values of x and Tmux for the modified models 
(except X for MMK) are comparable to those found earlier for the original models (§6 of 
BW). 



4.1. Comparing Models with Additional Statistical Tests 

Examining the 2-D K-S test results given in Table [5] we can say that the [P^D] 
and [z-D] planes can be reasonably fitted by the "varied" cases of the modified models. 
Six of these 9 planes not involving a (those three slices for each of the 3C, 6C and 7C 
comparisons) had 2-D probabilities > 0.1 for the MMK "varied" model; this is true for 4 V's 
in the MBRW cases. All of the 2-D P's of the "varied" MMK model are higher than those 
of the default MMK. When compared to the corresponding default versions, improvements 
are seen for 8 of 9 of the 2-D P's not involving a in the "varied" MBRW model. These 2-D 
results provide support for the superiority of the "varied" models (selected from 1-D K-S 
tests) in fitting the data. 

Comparing the "varied" cases of the two modified models themselves in the 9 planes not 
involving a, we see that 6 of the 2-D P's for the MMK model are higher than those of MBRW. 
Nonetheless, the a-related 2-D K-S probabilities are < 0.008 for every modified model. So, 
similarly to the original models, the modified models cannot fit any plane involving a. From 
the 2-D K-S probabilities we conclude that the MMK model is better (having the highest 
number of 2-D P's close to 1) in fitting the observational data than is the MBRW model. 
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Fig. 3. — The 3CRR, 6CE and 7CRS survey data, overplotted with the [P-D-z-a] planes 
for the 3C, 6C and 7C simulations of a good fit of the MBRW model. The initial ensemble 
(of size 4963343) is generated using x = 3.0, Tmux = 300 Myr; the power evolution is with 
the parameter variation of t^f = 100 yr, with the rest being their default values as in the 
MBRW model O. The 1-D K-S statistics for this case are in Tabled] (last but 6th entry). 

From the Spearman partial rank correlation analyses on the combined data of the 3 
surveys (Table [6]) we conclude that many matches to the data correlations are acceptable 
for the MBRW model, but they are less good for the MMK model. It is interesting to note 
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that the parameter variation cases which were the best fits (i.e., gave the highest combined 
probabihty, V[p^D,z,a]) from the 1-D K-S results, or the "varied" cases in Tables [5] and [6l are 
not necessarily the best fits according to the correlation analyses. For the MBRW and MMK 
models the default and the "varied" cases perform comparably, as 3 correlations are better 
in the default, and the remaining 3 are better in the "varied", models. 

Considering the signs of the four- variable coefficients of the combined surveys, the MMK 
model predicts a [-P-a] anti-correlation, a [D-z] correlation; its "varied" case also produces a 
[P-z] anti-correlation. All these are trends opposite to those seen in the survey data and to 
the other models, except that the MBRW model also predicts a [-P-a] anti- correlation. Both 
the MMK models are the only ones which produce the correct signs of the [-D-a] and [z-a] 
correlations of the combined survey data (the [-D-a] correlation is also exhibited weakly by 
the MBRW default case). However, this advantage is not meaningful as the MMK model 
gives very poor fits to the actual a-distributions. 

According to the correlation coefficient analyses involving the data from all the surveys 
(3C, 6C and 7C) the the MBRW model fits the data most closely, followed by MMK. In 
particular, the MBRW model provides the best fit to the key 4- variable correlation, rpp) ^a- 
This indicates that in the BRW model a growing hotspot is able to reproduce the P-D 
evolution (seen in 3C, 6C and 7C survey data) better than assuming a constant hotspot 
size (the original BRW model). Similar trends emerged when we examined the 3- variable 
correlation coefficients. 



4.2. [P-D-z-a] Planes 

We plotted planes through the [P-D-z-a] volume for the simulated surveys using the 
modified models, and compared them with the overall trends in the [P-D-z-a] slices of the 
observational data. The actual 3C, 6C and 7C data are plotted as crosses and the simulated 
data (3C, 6C and 7C virtual surveys) for the "best-fit" or "varied" parameter set for the 
MBRW model (those which give the highest total 1-D K-S probability) are overplotted as 
squares in Fig. [31 The Monte Carlo results for the "best-fit" parameter set for the MBRW 
model are shown in Fig. |H The main features of the [P-D-z-a] planes of the modified 
models are analogous to those of the original models in §6 of BW, so we discuss them only 
briefly, placing stress on any new features. 

In the [P-z] plane, all of our simulated surveys of all the modified models miss many 
of the low-P sources seen in the data. Too few low-2;/low-P sources are produced in all the 
simulated 7C surveys. There is underproduction of very high-2; sources {z > 2) in the 7C 



- 15 - 



3CRR Simulation 



6CE Simulation 



7CRS Simulation 




1 2 
Redshift, z 




1 2 
Redshift, z 




xn □ „ 

* *+ +.++ 



10 100 1000 

Projected Linear Size, D (kpc) 



10000 





29 ■ 




28 - 


'WHz 


27 - 


qT 


26 - 


o 
_I 


25 - 

24 



10 100 1000 

Projected Linear Size, D (kpc) 



10000 



_ 29 
J° 28 

N 

I" 

S 26 

i 25 

o 

_i 

24 



1 




10 100 1000 10000 

Projected Linear Size, D (kpc) 



Q. 10000 

Q 

a; 1000 

N 
</> 

TO 100 



#1* X 

A Aif^A □ □ 



1 2 

Redslnift, z 




1 2 
Redshift, z 



o 

Q. 


10000 ■ 






Q 






1000 r 


N 




CO 




near 


100 r 




_i 




■o 




10 - 


t) 






o 


1 ' 


qI 






A^ i^^^^ 



1 2 
Redshift, z 




1 2 
Redshift, z 




1 2 
Redshift, z 




1 2 
Redshift, z 



Fig. 4. — The [P-D-z-a] planes for the 3C, 6C and 7C simulations of a good fit of the 
MMK model. The initial ensemble (of size 3888492) is generated using x = 2.6, TMax = 150 
Myr; the power evolution is with parameter variation (5 = 1.6, the rest of the parameters 
having their default values of the MMK model (0). The symbols classify the sources into 
redshift bins as follows; Plus: < z < 0.5, Triangle: 0.5 < z < 1.0, Cross: 1.0 < z < 1.5, 
Square: 1.5 < z. The 1-D K-S statistics for this case are in Table [2] (9th entry). 

simulations and a similar, but less pronounced, trend is also present for 6C. 

From the trends in the [P-D] plane, the MBRW model overproduces large powerful 
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sources in 3C, and underproduces the large weaker sources. The 6C and 7C [P-D] planes of 
all the modified models show a closer match to the data, but the actual data is more scattered 
than the simulations. Given that any additional physics not included in the models would 
tend to broaden this distribution, this result is expected. 

There is P-D anti-correlation in the MMK model in all of the 3C, 6C and 7C simulations. 
Such P D evolution is also seen in the MBRW model where it is more pronounced in the 
6C and 7C simulations. An important improvement in the MBRW model is that the strong 
P-D anti-correlation apparent in even the "best-fit" BRW model is diluted when modified 
to incorporate a growing hotspot. 

In the [D-z] plane, the MBRW model still overproduces the larger 3C sources. The 
MMK model's "best-fit" case [D-z] planes seem to be a good fit to the data (especially 
the 3C [D-z]), with about the right amount of D-z anti-correlation. The MBRW model 
produces a weaker D-z anti-correlation than does the data in 3C, but one stronger than 
in the data in 6C and 7C. An explanation of the D-z evolution as an effect of the "youth- 
redshift degeneracy" was given while discussing the original model planes in §6 of BW. From 
these trends we conclude that the MMK model is the best fit (by eye) to the [D-z] planes 
of the 3C, 6C and 7C data. 

All the other characteristics examined in the simulations show tight correlations with 
the spectral index. The data shows much greater scatter in a than is seen in any of the 
models. The spectral index distributions of all the model simulations are similar to those in 
BW. Since the a statistical fits are very poor, we do not discuss them any further. 



4.3. Alternative RLF 



The iGrimes et al.l ( 120041 ) RLF is tested using different values of x and TMax, but with 
the default values of radio lobe power evolution model par ameters of the different models. 
When comparing th e simulations using iGrimes et al.l (120041 ) RLF with respect to those with 
Willott et al.l (120011 ) 's (with the same model parame ters otherwise), most of the models give 



better 1-D K-S statistics using IWillott et al. 



(j200l|)'s RLF. The KDA and MK models are 



better for 5 out of 6 comparisons, the MBRW model is better for 4 out of 6 total. The MMK 
model perfor ms equally us i ng bo th RLFs. The BRW is the only model which give better 
results using IGrimes et al.l (120041 ) RLF for all 5 comparisons we made. 



S o we conclud e that the KDA, MK and MBRW models give better fit to the data 



using IWillott et al.l (l200ll )'s RLF. From the runs we performed the MMK m odel produces 
comparable results with both the RLFs. The BRW model is better fit with the IGrimes et al. 
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( 120041 ) RLF, but in an absolute comparison with the other models (KDA, MK and MBRW), 



the fits are still poor. 

All the spectral index fits are still very poor using the alternative RLF. So we can say 
that altering the redshift birth function has little effect on the major drawback of these 
models: the mismatch of the spectral index behavior between data and simulations. 



5. Relevant Volume Filling Fraction 

The volume of our "Relevant Universe" is taken to be the volume of the cosmic baryons 



1999. 


2006; 


Dave et al. 


2001) 



< T < 10^ K (e.g., ICen fc Ostriker 



large fraction of the baryons (~ 40 — 50% by mass) in the universe at the present epoch. The 
WHIM permeates the universe as extended large-scale filamentary structures, the junctions 
of which are the sites of galaxy and cluster formation. 

To calculate the volume fraction of the relevant universe filled by RGs we first consider 
that the initial ensemble of sources generated to yield the simulated surveys provides a good 
estimate of the actual population of RGs at different redshifts since, as discussed above, 
those Monte Carlo simulations do produce decent fits to the P-D-z distributions. Among 
the millions of RGs generated in such an ensemble, only a few tens to hundreds are detected 
in the simulated surveys. The key point is that this is also what happens in reality: among 
the millions of RGs born through the quasar era, we can now detect only a few in our fiux- 
limited and redshift-complete radio surveys. Severe energy losses (adiabatic, synchrotron 
and IC losses) conspire to force the vast majority of all RGs to fall below our observational 



flux limits, as discussed in detail in earlier work (iGopal-Krishna. Wiita. &: Saripallil Il989 



Blundell fc Rawlind Il999l : BRW; GKWOl; GKW03a,b; GKWO; GKWB04). Such power 



losses are evident in the simulations from the steeply falling [P-D] tracks of §3.1[ 

Next, the universe is divided into redshift bins (shells), and the relevant volume fraction 
is calculated in each 2;-bin. Let the minimum and maximum redshifts of a shell be Zmin and 
Zmax- The bin-width is taken as /S.z = 0.02 = {zmax — Zmin)- The mid-redshift of a bin, 
Zmid = {zmin + z^nax) /2, is cousidcrcd as the epoch of that bin. Then, the distribution of 
the volume fraction is estimated as a function of redshift, which is then integrated over the 
entire quasar era epoch to get the total volume contribution of several generations of RGs 
in the universe. 
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5.1. Volume of WHIM in the Universe 



We use a consensus flat, dark-energy dominated universe (ISpergel et al.l l2003l . l2006l ). 
witli Hubble constant Ho = 71 km s~^ Mpc~^, matter density parameter Qm = 0.3, and 
vacuum energy density parameter Q\ = 0.7. The como ying volume over all-sky {An sr) in a 
redshift shell of the universe between z^in and z^ax is ( HogglEoOol ). 

3 

(7) 



^^=-f (i) ^-^ 



where, 



X2,l 



[nMii + zf + nA] ^'^dz. 



Our big ensemble (initial population) explicitly detects sources over the 3CRR sky 
survey area, Areas^ = 4.23 sr. So the RG population from which these sources are detected 
lies within a smaller comoving volume extending over only the 3CRR sky area, since in these 
simulations the number of sources can be taken as proportional to the sky area over which 
they are surveyed. If a simulation detects A^sim(3C) sources, where there are Ai'samp(3C) = 145 
sources in the real 3CRR survey, then the 3C detection ratio is written as, Ratiosc = 
Nsim(3C)/Nsamp{3C)- Heuce the comoving volume over the effective 3C survey area is, 

Areasc x Ratiosc 



co,3C(eff) 



(9) 



The effective comoving volume of the 2;-shell is then converted to the proper volume it had at 
that epoch, as AV^^op = A Vco,3c(efr) / (1 + Zmid)'^ ■ The effective relevant volume of the 2;-shell 
is then the fraction of the proper cosmological volume of the shell occupied by WHIM. 

The WHIM volume fr action, is adopted fr om the large-scale cosmological simulations 
of ICen &: Ostrikerl (119991 ). ICen fc Ostrikerl (120061 ) give an improved WHIM fraction calcula- 
tion by explicitly including galactic superwind feedback processes, but for our purposes there 
is no significant difference from their previous results. Hence we obtain the final "relevant 
volume of the universe" inside a z-shell, AVwhim = ^^Vprop- 



5.2. Radio Galaxy Volumes and Relevant Fraction 



We calculate RG volumes by assuming that the RGs are cylindrical in shape with total 
length D{t) (Eq. 4 of BW), at an age t. The axial ratio, Rt gives the ratio of the source 
length and its width (or diameter). The volume occupied by a RG at an age t is, 

2 _7^/.N3 



VRcit) 



TC 



2Rt 



71 



ARl 



(10) 
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To get a conservative estimate, in all the volume computations Rt = 5 is used, irrespective 
of the model (unless otherwise noted). This value appears to be an upper bound to the 
average axial ratio based on observations (GKWOl and references therein). The difference 
between this choice of Rt = 5 used to calculate the volumes, and that in the KDA model, 
where Rt = 1.3 gave the best fit to the [P-D-z] planes, is noteworthy. If later work shows 
that Rt < 5 is indeed preferable, then a typical radio galaxy volume Vrg(^) will be bigger, 
by up to a factor of ~ 15, thus more strongly favoring the picture of substantial cosmological 
impact of RGs. 

Following the conservative bent of GKWOl, we only consider the contributions of the 
more powerful FR II RGs despite the fact that the weaker FR I RGs are much more numerous 
in our local universe and some of them are also seen to extend for hundreds of kpc. There are 
several good reason s for making this ap proximation, although we cannot precisely quantify 



its effect. First, the lWillott et al.l (120011 ) RLF indicates that the numbers of FR I and FR II 
sources were much more comparable during the quasar era on which we are focusing, with 
the lower luminosity sources hardly evolving for z > 0.7 while the higher luminosity sources 
continue to rise in density with z until z > 2. Second, it is very likely that the efficiency 
with which je t thrust is converted into rad i o flux in FR I's is substantial ly higher than in 



FR IPs fe.g.. iGopal-Krishna fc Wiitalll99ll : iBaum. Zirbel. fc 0'Dealll995[ ). so that a jet of 



equal power would make a brighter FR I source, thereby skewing the RLF of FR I's upward. 
But since the typical FR I has much lower radio flux it is being energized by very much 
weaker jets than those powering a typical FR II. Hence, on average, the volumes enclosed by 
individual FR I's will be substantially lower than those of the FR II RGs considered here. 
In addition, since the expansion of FR I's is less rapid they produce weaker (if any) shocks 
for shorter times and will have less ability to trigger star formation. 

Let the cosmic times of the boundaries of the z-shell be denoted as tin and tout (cor- 
responding, respectively, to Zmin, Zmax)- We assume that all the RGs in an ensemble hve 
out to their full lifetime, TMax- So, while the vast majority of FR II radio sources ever born 
are too faint to be detected now, they do expand as long as the AGN is feeding the jets 
and lobes and hence contribute to filling part of the universe. We count all the RGs in the 
simulation initial ensemble which have any portion of their lives falling in the time range of 
that z-bin. The volume contribution of all the RGs which are intercepted by a z-shell are 
then added to get the total RG volume, AVrg- The relevant volume fraction in a z-shell is 
/S.l{z) = AVrg/^Hvhim • Integrating At{z) over z gives this volume fraction as, 

^early 

AL{z)dz, (11) 



where Zeariy is the earliest redshift of a source in the initial ensemble. 
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Fig. 5. — Relevant volumes in the Universe: WHIM volume over the effective 3CRR sky 
area, and total RG volume, for the model with parameters noted as a function of redshift. 



There is another significant factor which must taken into account to estimate the total 
fraction correctly; this arises from the contributions of the several generations of RGs during 
the quasar era (QE). The total volume filled by the multiple generations of RGs in the 
universe over the whole QE was roughly taken into account in GKWOl. They considered 
the length of the QE as tqE ~ 2 Gyr, and the maximum age of radio sources Tjvfax = 500 
Myr. We also adopt this rather large maximum length of AGN activity in some of our 
simulations, as it was suggested by BRW. Ages of that order are supported by a variety of 
recent observational st udies mentioned in §2 of BW; the most recent modeling of RG ages by 



Machalski et al.l (120061 ) also implies high values for Tuax- GKWOl argued that every place 
in the universe could have been potentially affected by tqE/TMax = 4 generations of RGs 
during the entire QE. So they multiplied the mean of the corrected RLE by (tgs/^A/aa) to 
get the total proper density of intrinsically powerful radio sources in the universe. 

In our simulations we obtain the total fraction by adding the values of /S.l{z) several 
times in intervals of TMax over the entire QE. The length of the QE is obtained from the 
temporal length of the epoch for which Al{z) > 5% of its peak value. Starting from the 
high-2; end-point of the QE, values of Al{z) are computed at intervals of Tmqx and summed, 
until the \ow-z end-point of the QE is reached or exceeded. This addition is done several 
times; each time the starting point is chosen differently by going back or forward from the 
original starting point by integral multiples of 50 Myr. The summed Al{z) obtained from 
these several additions (each starting from a different cosmic time), are then averaged to get 
the mean total "relevant volume fraction", (, of the universe filled by generations of radio 
galaxies during the quasar era. 
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Fig. 6. — Relevant volume fractions of the universe filled by radio galaxies from several 
model simulations. The vertical alignment of the labels correspond to those of the model 
plots. 



The relevant volume fraction, C, was computed for a subset of the model simulations 
done with the 3 main models, KDA, BRW and MK, as well as for the modifications we have 
considered in the present work. 

We now consider the model parameters which determine the relevant volume frac- 
tion. The distribution functions of RG z (Eq. [6j zo,az), and of jet power (Eq. 2 of BW: 
X, Qmin, Qmax), accordiug to which an initial ensemble of sources are generated following the 
presciption from BRW, along with the maximum age, Tmqx, are the parameters which are 
model-independent in the sense that they do not depend on the RG lobe power evolution 
models. The RG volume, Vrg(^) (Eq- [ID]) depends on the models through the linear size, 
D{t) (Eq. 4 of BW), which explicitly involves the ambient density parameters (Eq. 3 of BW: 
Po,ao,j3). The other model-dependent factor is the detection ratio ( §5.ip which is obtained 
from the number of sources in the simulated surveys. 



5.3. 



Results and Discussion 



5.3.1. Cosmological Volumes vs. Redshift 



In Fig. the dashed curve is the final relevant WHIM volume AVwhim ( §5.11) . and 
the solid curve is the total RG volume (AVrg, §5.21) . for a case of the KDA model. The 
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RLF or redshift distribution from which the sources in the initial ensemble are drawn is a 
gaussian, so the maximum number of sources in the initial ensemble are born near the peak 
at zq = 2.2. From Fig. Owe can see that the total RG volume, AVrg, peaks at z ^ 2 < zq. 
This is because the majority of sources born at zq remain active for T^ax = 500 Myr, and 
are thus counted in several later 2;-bins. They contribute to the RG volume in increasing 
amounts as they grow in age, until the maximum age, Tmux, after which they are assumed 
to no longer do so. Their combined increasing contributions at later times make AVjig peak 
at a z < Zq. This peak redshift of AVrg should be around the cosmic epoch corresponding 
to t{zo) + Tmux, as that is when the largest number of sources in the population reach their 
maximum volumes. 

On the other hand. Fig. [H] shows that the relevant volume fraction, Al{z), peaks at 
z > Zq. The distribution of WHIM volume, AVwhim (Fig. E]), can be invoked to explain this 
result. We see that AVwhim rises sharply from z ~ 3, until it reaches a peak at z ~ 0.7, 
because of the trends of proper volume of z-shells in the consensus cosmology. When the 
ratio of AVrg to AVwhim is taken to get Al{z) at 2; ~ zq, AVrg is divided by a volume 
AVwhim which decreases with increasing redshift. 

A representation of how the volume contribution of multiple RG generations are added 
to get the total cumulative fraction over the entire QE is given in Fig. [3 The solid black 
curve in the figure is the volume fraction Al( clS cl function of redshift. The symbols (of a 
single type) plotted on it are the values of Ai{z) which are picked at intervals of TMax = 500 
Myr over the QE, and added. The different plotting symbols denote the different starting 
points for the added fractions, which were finally averaged to get the cumulative fraction, 
^. For this model of Fig. [7] (BRW default using initial ensemble with x = 2.6, Tmux = 500 
Myr), the quasar era spans the redshift range zqe = 3.52 — 1.16, or the cosmic time range 
1.74 — 5.10 Gyr, corresponding to a quasar era of duration Iqe = 3.36 Gyr. Hence there are 
contributions from ~ 7 generations of RGs in the case where Tmux = 500 Myr. The final 
relevant fraction results for this model are: 

L = 0.0123, C = 0.0301. (12) 



5.3.2. Comparison Between Models 

Fig. E] shows plots of the relevant volume fraction, Al{z), as a function of redshift, for 
some of the model simulations. Table [7] gives the relevant volume fraction results for the 
models. The results for a particular parameter variation of a model are given in each row. 
Most of the column headings were defined earlier. Column 6 gives A^, the normalizing factor 
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Fig. 7. — Total relevant volume fraction of the universe filled by RGs by adding the volume 
contributions of multiple generations of RGs over the quasar era. This is for the BRW 
simulation with default model parameters, for the initial ensemble of size 1561417 generated 
using X = 2.6, Tmux = 500 Myr. 

in the redshift distribution (Eq. [6]) which was used to generate the initial population; 
is the factor by which [Vc{z = 0) x p{zhirth)] (paragraph 5, §2 of BW) is multiplied to get 
the number, Nbom, of radio sources born within the relevant comoving volume, Vc- From 
Table [7] we see that the fractions l and ( vary significantly both between the models and 
for different parameter sets within the same model. Quite a wide range of relevant volume 
filling factors can be produced by the considered ranges in model parameters. 

For the "default" case, RGs in the BRW simulations cumulatively fill ~3% of the relevant 
universe, and ~3.4% of the relevant volume is filled in the MBRW model. This number is 
~4.5% for the KDA, ~7% for the MK and ~13.5% for the MMK models. At the same x 
and Tmox, and using default model parameters, the MK model gives the largest relevant 
fractions, followed by the KDA and finally by the BRW models. This is because a larger 
initial ensemble needs to be generated in the MK model to yield the same number of detected 
sources compared to the other models. 

Still, the runs involving parameter variations corresponding to "1-D K-S best-fit" s of the 
BRW and KDA models give higher fractions (by ~ 4 — 5 times) than do the corresponding 
"default" models, for the same x and Tmux- This is because the BRW best-fit is with ao = 7.5 
kpc < ao (Default), and KDA best-fit is with po = Po (Default) /2, both of which have the effect 
of increasing D{t) and thence the RG volumes. However, when the favored lower values of 
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TMax are employed, the cumulatively filled fractions are 0.015 < C < 0.07. The relevant 
volume fractions for these most preferred (with respect to K-S statistics) parameter sets are 
all very low compared to the estimate in GKWOl, ( = 53%. Still, for certain less likely 
parameter values this fraction goes as high as 16% for BRW, 55% for KDA, and 20% for 
MK model. 

The relevant fraction is greater for higher Tmux (for same x), as is evident from the 
expression for the RG volume, which scales as For the same T^ax, the fraction 

is higher with x = 3.0 than with x = 2.6. This might seem counter-intuitive, as with a 
steeper jet power distribution the sum of volumes occupied by the same number of sources 
should be smaller. The mean RG volume at maximum age, {V{TMax)) is indeed smaller (by 
a factor of 1.21/1.44 = 0.84), as discussed later ( §5.4p . Nonetheless, a larger volume fraction 
at higher x can be explained as many more sources must be generated using the steeper 
slope {x = 3) to yield numbers of sources in the simulations comparable to those in the real 
surveys. Explicitly, to get the same Ratiosc, larger ensemble sizes (by 1.5 — 3.5 times) are 
required for x = 3 than for x = 2.6, and the increase in ensemble size more than offsets the 
smaller mean volume. 



5.4. Comparison of Results with a Previous Estimate 



Gopal-Krishna fc Wiital (120011 ) performed a preliminary calculation to find the relevant 
volume fraction added over several generations of radio sources during the quasar era. From 
the results in §5.31 we see that the relevant fraction obtained in this work is considerably 
smaller than the fraction estimated by GKWOl. For the BRW default model (simulation 
done with initial ensemble of size 1561417 generated using x = 2.6, Tmux = 500 Myr) we 
obtained ( ~ 0.03 (2nd entry of Table [7]), whereas using published graphs in BRW for the 
same model GKWOl obtained ( ~ 0.5. 

There are six differences between the calculations, which taken together, can account 
for most of the discrepancy between the two calculations. 



[1) We adopted the newer consensus cosmology (e.g. ISpergel et al.ll2006l ). GKWOl used 



cosmologies with Hq = 50 km s ^ Mpc ^, and with either Qm = or 1, but with = 0. 

(2) In the model simulations performed in BW and this work, we used Qmin = 5 x 10^^ 
W as the minimum jet power, following BRW (as described in §2 of BW). GKWOl took 
the effective lower limit of Qo as Qm = 7.5 x 10^''' W (their §2.2), which they inferred (by 
observing the BRW [P-D] tracks in their Figs. 13 and 14) to be the minimum power a 
source must have in order to appear in the BRW data set. A higher minimum Qo means 
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that the RG jets will be, on average, more powerful, thus making the total RG-volume and 
the relevant fraction higher. 

(3) We obtained the average radio galaxy volume at maximum age Tmux = 500 Myr to 
be 1.44 Mpc^ whereas GKWOl obtained a value of 2.1 Mpc^ (their §2.3), using the same 
BRW default model parameters (except for the different Qmin choice). Clearly, {V{Tniax)) 
should scale as Q^in ^\ since Qmax ^ Qmin- So GKWOl's {V{TMax)) should be higher than 
ours by (7.5/5)^/^^"^^ = 1.42. 

(4) To get the total proper density of RGs, GKWOl multiplied the peak of the corrected 
RLF by [tQE/TMax) = 4 generations of RGs in their §2.2. We consider the contribution from 
multiple generations of RGs by adding the volume fractions Al{z) in intervals of TMax over 
the entire QE ( §5.21) . so we more precisely take into account the distribution of Al{z) vs. z. 
As can be inferred from Fig. [71 simply multiplying the peak fraction (highest Al{z) in Fig. [7]) 
by 4, which was done in GKWOl, gives an overestimate of C by a factor of 0.038/0.030 = 1.27. 

(5) We interpolate d the fractional WHIM volume in the universe, which varies with z, 



from the calculations of ICen fc Ostrikerl (119991 ). This WHIM volume fraction decreases with 



increasing redshift, starting from ~ 0.095 at the present epoch, to reach ~ 0.01 at z ~ 3. 
GKWOl considered a contant WHIM volume fraction of 0.03 at all redshifts of the QE. 
For our computation, the higher WHIM fraction at low- 2; dominates over the lower WHIM 
fraction at high-z. 

(6) In converting from comoving to proper volumes GKWOl used a value of z = 2.5. 
We integrated over each value of z and had an effective average value of z ~ 2.2. This 
value is basically the peak, zq, of the Gaussian redshift distribution (Eq. [6]) of the radio 
sources. This difference causes the GKWOl result for ( to exceed ours by another factor of 
~ (1 + 2.5)7(1 + 2.2)3 = 1.31. 



6. Conclusions and Future Work 

We have performed comprehensive quantitative tests of extensively explored modifica- 
tions to the BRW and MK models for FR II RG evolution which allowed the sources' hotspot 
sizes to grow with age. We often found adequate fits to the [P-D-z] distributions for each 
model for each of the Cambridge catalog subsamples (3CRR, 6CE, 7CRS). But we cannot 
locate any parameter sets which provide good simultaneous fits to all three catalogs and to 
all four of these observables [P-D-z-a]. Of particular concern are the spectral indices, where 
none of the models provides an adequate fit. We note that all of these models only produce 
a single gross spectral index and that future models must improve upon the treatment of 
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the physics that leads to average emissions at different frequencies. 

Prom the power vs. hnear-size tracks we see that the BRW [P-D] tracks are the steepest 
among all six models, while the MBRW tracks are significantly less steep. This causes the 
most significant change found in the model performances after modification, in that the 
MBRW model gives substantially better K-S statistical fits to the data than the BRW model. 

All the models produced better statistical fits with the slope of the jet power distribution 
set to X — 3, except the MMK model which performed better with the default value from 
BRW, X — 2.6. Considering the active lifetime of the AGN for which the jets feed the lobes, 
we found that all the models gave better fits with TMax between 150 - 300 Myr; the high 
default value of 500 Myr is disfavored. 

The KDA, MK, MBRW and MMK models all perform comparably in terms of producing 
high values of total 1-D K-S probabilities. Prom the 2-D K-S test results the "varied" cases 
of most of the models can produce adequate fits to the [-P-^], [-P^-D] and [z-D] planes. 
Sorting all the models in descending order of the number of non-a 2-D P's greater than 
any other model we have: KDA, MK, MMK, MBRW, BRW. From the 4-variable Spearman 
partial rank correlation coefficient analyses, wc find that the original KDA model can match 
the survey data correlations very closely (at least for P, D and z), followed by BRW and 
MBRW, then finally the MK and MMK models. 

The existing redshift complete radio source catalogs can limit the allowed parameters for 
each of the models individually and can show that the K2000 model is not a good fit to the 
ensemble of RGs. However, they are too small to make a clear distinction among the models 
we have tested. A significant range of the models can provide adequate descriptions of the 
[P-D-z] distributions. Only by obtaining redshift-complete samples covering significantly 
larger portions of the sky than those studied in the 6C and 7C surveys would the data have 
enough power to positively discriminate among the models. 

A major goal of this work has been to calculate what fraction of the "relevant universe" , 
or the large-scale filament-structured WHIM volume, do the PR H RGs born over the quasar 
era cumulatively occupy. We found that quite a wide range of relevant volume filling factors 
can be produced, with certain choices of parameters in our simulations producing relevant 
filling factors as high as ~ 20% and even ~ 50%; however, none of these models provide 
good K-S statistical fits to the observations and hence these rather high filling factors are 
unlikely to be realistic. For both the default and 1-D K-S "best-fit" cases of the models 
using the best values of TMax the cumulative volume filling fraction of several generations of 
RGs varied within the range ~ 1.5 — 7%. These volume fractions are significantly smaller 
than the preliminary estimate of GKWOl. 
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We infer that FR II RGs probably cover ~ 5% of the WHIM volume cumulatively over 
the quasar era. Such filling factors were produced by models providing the best statistical 
fits to the observations, which are quite good fits to the P-D-z distributions, so they are 
likely to be close to the most probable "real" fractions. Still, despite the large number of 
models we have examined they are not numerous enough to allow a proper error estimate on 
this key quantity. The fractional "internal" errors on these filling factor values, as obtained 
by computing results for several of the same models and same parameter values but with 
different initial random seeds, were in the range 0.017 < cr/{C) < 0.10. 

We conclude that the expanding radio galaxies born during the quasar era can play a 
significant, though probably not dominant, role in the cosmological history of the universe 
through the triggering of extensive star formation and the spreading of magnetic fields. 



Recently iBlundell et al.l (120061 ) provided observational evidence for the discovery of low- 
energy cutoff of particle acceleration in the lobe of a giant FR II RG. They obtained a value of 
lmin{hs) ~ 10'^ as the minimum Lorentz factor of particles in the hotspot, substantially higher 
than the values in the models where (1 — 10). Investigating models using such (tentatively) 



observationally supported higher 7, 



min{hs) 



values will be a worthwhile venture. 



In a recent study iKawakatu fc Kind (120061 ) described the dynamical evolution of the 
hotspots of radio loud AGN using a model significantly different from that of KDA. Incor- 
porating such evolving hotspots into the models of radio lobe power evolution we analysed 
would be an interesting modification to the models worth exploring. 

A potential indicator which can provide an excellent test of whether RGs do really 
trigger galaxy formation, is the 3-poi nt correlation functio n between radio galaxies in large 
scale galaxy (redshift) surveys (e.g., iBorderia et al.l 1 199 if ). If our RG impact scenario is 
robust, then there should be a bias in this correlation function along the direction of radio 
lobes of the RGs, because more galaxies are formed along the radio-axis triggered by jet/lobe 
expansion, as compared to directions perpendicular to the radio-jet. 

The models investigated in this work and in BW predict the power from the radio 
lobes only. A natural extension involves the question of whether the same models also fit 
deeper radio catalogs if we take into account the relativistically Doppler boosted core/jet 
emission. By incorporating the beamed core emission, investigations of simulations of large 
scale radio surveys containing many thousands of sources can be done . Some such deep 
surveys are Faint Images of the Radio Sky at Twenty-cm, FIRST ( iBecker et al.l Il995l ) , the 
Westerbork Nor thern Sky Survey, W ENSS jRengelink et al.lll997h and the NRAO VLA Sky 



Survey, NVSS (j Condon et al 



9981). which can be made adequately complete in redshift 



through optical identifications (llvezic et al.ll2004l ) from the Sloan Digital Sky Survey, SDSS 



( lYork et al.l |2000| ). If simulations can be performed to predict thousands of sources, the 
possibihty of successfully incorporating a multi- dimensional statistical test becomes much 
greater. 
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Tablc 1. MBRW Model: 1-D KS Statistics for Selected Parameter Variations 
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Table 1 — Continued 
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ouu 




\J\J 


794 


0'i82 


u.uuu 


9 1 1 p-09 

. X X C U tJ 


2 26 






7C 


2.36e-04 


0.741 


1.35e-04 


1.46e-05 




3 n 

o . w 


f 

/In /lo 


3C 


022Q 


036 


061 6 

W . WU X u 


3 22p-1 1 

O. ^^\^ X X 


656 


3nn 

OUU 


4Qfi3343 

tiC/UOOtiO 




41 




287 




783 

U. / OO 








A AO ^ 


n 1 oo 


z.ooe-Oo 


r? Ad^ Cl^ 




3.0 


Tc = 5/3 


3C 


0.00283 


0.00201 


0.0105 


0.0361 


1.13 s 


300 


4963343 


6C 


0.987 


0.00293 


0.764 


0.0261 


1.14 s 






7C 


0.698 


0.00197 


0.569 


0.00528 




3.0 


^min{hs) — 10 


3C 


0.253 


1.36e-04 


0.579 


6.02e-ll 


2.31 


300 


4963343 


6C 


0.831 


0.227 


0.716 


7.62e-08 


2.81 






7C 


7.70e-05 


0.584 


1.35e-04 


0.00117 
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Table 1 — Continued 



X 

Tuax^ 


Model 
Ensemble Size 


Survey 


V{P) 


V{D) 


V{z) 


V{a) 




3 n 

Kj .yj 


^ /I, N — 10^° 

imax(ns) 


3C 


177 

U . i 1 i 


^ . I U\j uo 


171 

U. X 1 i 


3 29p-1 3 


2 29 


ouu 




fin 


U.t^^O 


U.ZiO'4: 






ni 

o.ux 






7C 


0.00115 


0.932 


0.00111 


4.56e-04 




3 n 


ry inl6 

Imax(hs) 


3r; 


n 227 


(J . o 1 wo 


n 22(1 


5 35p-14 

O . OO t> Xt: 


1 nq 


3nn 


4Qfi3343 

tiC/UOOtiO 


U 


n ^33 

u . ooo 


n DDI 58 

U .UU X oo 


f) 272 

U. 1 




1 21 

X . Z X 






7C 


0.00455 


0.199 


0.00287 


0.00468 




3 n 


^ — Zi.UU± 


3C 

OV7 


OOSQS 

U.UUOc/O 


1 5Qp-05 

± .Ot7C UO 


01Q1 


3 QSp-16 

o.zjokj ^-\J 


31 5 
u.oxo 


3nn 


4Qfi3343 

tic/UOOtiO 




Ofifil 

U.UUUi 


001 58 

\j .yjyj A-OO 


00675 
u.uuu / o 


2 33p-07 

Zi .OOC U 1 


490 






7C 


0.0353 


0.286 


0.0344 


0.0198 




3 n 

o.u 




3r! 

Ov^ 


n 253 


2 3fip-n4 

.OUC Ut: 


451 

U . 'i O -L 


1 20p-1 1 

X . Z U C X X 


1 7Q 

X . 1 C/ 


ouu 


^t/UOO^O 


Uvj 


n 794 

U. 1 z,^ 


U .UOU 1 


u. oou 


7 9zLp HQ 


9 OR 
z.uu 






7C 


2.34e-05 


0.412 


7.56e-05 


0.00534 




o.u 


^ _ O QQQ 
JJ — Zj. cfcfcf 


3r! 

Ow 


U . J- c/O 


J- .rt J-C Ut: 


n 371 

U. O 1 X 


5 OQp-12 

O. Uc/C xz 


1 7Q 


3nn 

ouu 


4Qfi3343 


UV7 


SfiQ 


020fi 

U.UZiUU 


1Q4 

U. ±C/t: 


7 24p-0Q 

/ .Z'4:C Uc/ 


2 34 

Z . Ort 






7C 


8.02e-05 


0.879 


1.35e-04 


0.00468 




3 n 

o.u 




3r! 


n 1 30 

u . xou 


2 7np-n5 


0532 
u .uooz 


2 83p-1 5 

Z. OO" xo 


743 

U . 1 t:0 


3nn 

ouu 


4Qfi3343 

tiC/UOOtiO 




50Q 


n nnfi22 


n lOQ 


7 fi2p-08 


8fi3 
u.ouo 








0.03/ i 


o.iyo 


O.Oibi 


n ni no 
O.OiyZ 




3.0 


tbs = 10' yr 


3C 


7.23e-04 


1.90e-07 


6.75e-04 


1.28e-07 


0.397 


300 


4963343 


6C 


0.389 


0.0110 


0.0667 


2.11e-09 


0.477 






7C 


0.0113 


0.121 


0.0166 


1.46e-05 




3.0 


hf = 0.01 yr 


3C 


0.0978 


1.02e-07 


0.130 


2.79e-15 


0.809 


300 


4963343 


6C 


0.641 


0.00150 


0.188 


7.62C-08 


0.857 






7C 


0.00481 


0.0770 


0.00458 


0.00500 
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Table 1 — Continued 



X 


Model 












'P[P,D,z,a] 


rp a 
■J- Max 


Ensemble Size 


Survey 


ViP) 


V{D) 


V{z) 




V[P,2D,z,a] 


3.0 


hf = 10 yr 


3C 


0.290 


8.81e-06 


0.280 


2.91C-15 


2.06 


300 


4963343 


6C 


0.771 


0.0609 


0.945 


7.62e-08 


2.47 






7C 


4.03e-04 


0.603 


6.70e-04 


0.00976 




3.0 


hf = 100 yr 


3C 


0.253 


8.37e-05 


0.622 


8.34e-15 


2.66 


300 


4963343 


6C 


0.869 


0.234 


0.862 


7.62e-08 


3.34 






7C 


6.75e-04 


0.879 


3.99e-04 


0.00481 




3.0 


hf = 10=^ yr 


3C 


0.198 


2.36e-04 


0.579 


6.44^11 


2.05 


300 


4963343 


6C 


0.878 


0.00300 


0.862 


7.62C-08 


2.23 






7C 


3.96e-04 


0.286 


1.35e-04 


0.00500 




3.0 


hf = 100 yr 


3C 


0.113 


0.00379 


0.483 


2.80e-ll 


2.11 


300 


lmin(hs) = 10 


6C 


0.869 


0.159 


0.716 


7.62e-08 


2.62 




4963343 


7C 


3.52e-06 


0.668 


4.28e-05 


0.00917 




3 n 


f,j. — inn vr 


3C 


n n624 


2 66p-n6 

^ .yjyjK^ yjyj 


n 247 


1 76p-n6 


1 58 


300 


Oq = 15 kpc 


6C 


0.416 


0.0992 


0.556 


2.11e-09 


2.23 




4963343 


7C 


0.00115 


0.971 


7.67e-05 


5.49e-04 




3.0 


'ymin(hs) 10 


3C 


0.115 


3.97C-07 


0.192 


3.54e-09 


1.17 


300 


Gq = 15 kpc 


6C 


0.564 


2.50e-05 


0.408 


7.24e-09 


1.42 




4963343 


7C 


6.84e-06 


0.412 


7.67e-05 


3.76e-05 




3.0 


hf = 100 yr 


3C 


0.147 


1.07e-07 


0.192 


0.0385 


1.17 


300 


KDA Env. '=(/3,ao,Po) 


6C 


0.684 


0.0626 


0.556 


6.89e-07 


1.22 




4963343 


7C 


3.45e-06 


0.0163 


1.27e-05 


0.0103 
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^TMax in units of Myr. 



^All other parameters are as in t he MBRW model follo wing iBlundell et al.l (119991 ). with 
hotspot size increasing according to lJeyakumar &: Saikial (l2000l ). 

■^Parameters defining the external environment density profile are set to those of the KDA model: 
p = 1.9, ao = 2 kpc, po = 7.2 x 10-^2 kg m'^. 

"pi = Po (Dcfauit)/2 = 8.35 X 10-24 kg m-3. 

V2 = 2 X Po (Default) = 3.34 X 10-23 kg m-3. 

V3 = 4 X Po (Default) = 6.68 X lO^^s kg m'^. 

^Numbers of sources detected in some of the simulated surveys are considerably smaller than 
in the real surveys, so the 1-D K-S statistic does not hold much significance. 
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Tablc 2. MMK Model: 1-D KS Statistics for Selected Parameter Variations 



X 


Model 












'P[P,D,z,a] 




Ensemble Size 


Survey 


V{P) 


V{D) 


V{z) 


Via) 


V[P,2D,z,a] 


2.6 


Default ^ 


3C 


0.0931 


0.00650 


0.124 


2.80e-45 


0.727 


500 


4397469 


6C 


0.105 


0.0283 


0.575 


1.82e-21 


0.767 ^ 








0'^56 


O'^OQ 


T^O 

U. -LOU 


7 q9p_99 




2.6 


Default 


3C 


0.0491 


0.521 


0.535 





2.12 


150 


3888492 


6C 


0.624 


0.102 


0.506 


1.83e-24 


2.88 








00187 


"^01 


04Q1 


6 1 '^p-91 




2.6 


Default 


3C 


9.76e-04 


0.216 


0.421 





1.64 


250 


4195764 


6C 


0.0935 


3.44e-04 


0.681 


1.83e-24 


2.36 








Z/ .OOt: u<J 


n 897 


n ni 1 1^ 






2.6 


Default 


3C 


0.00152 


0.0125 


0.421 





1.72 


300 


4342468 


6C 


0.641 


0.252 


0.0709 


1.83e-24 


2.34 








0105 


739 


935 


6 1 3p-91 




3.0 


Default 


3C 


1.21e-07 


0.0366 


1.47^08 





0.558 


150 


4861474 


6C 


0.212 


0.371 


0.0265 


2.99e-26 


0.911 








00'^41 


1 '\A 


0541 


7 Q9p-99 




3.0 


Default 


3C 


2.30e-04 


0.195 


6.52e-05 


1.58e-31 


0.410 


500 


4886474 


6C 


0.154 


0.0125 


0.00722 


8.88e-16 


0.611 






7C 


0.145 


0.128 


0.0532 


6.13e-21 




2.6 


KDA Env. 


3C 


0.00553 


0.00253 


0.485 





1.55 


150 


(3, ao, po 


6C 


0.157 


0.485 


0.556 


1.83e-24 


2.15 




3888492 


7C 


1.24e-05 


0.483 


0.00180 


2.79e-19 




2.6 


f3 = 1.0 


3C 


0.00368 


8.37e-05 


0.474 





1.27 


150 


3888492 


6C 


0.331 


7.50e-04 


0.498 


1.83e-24 


1.54 






7C 


4.26e-05 


0.434 


0.00710 


8.39e-22 
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Table 2 — Continued 



X 


Model 












T^[P,D,z,a] 




Ensemble Size 


Survey 


V{P) 


V{D) 




V{a) 


V[P,2D,z,a] 




/9 — 1 6 




n 1 21 


n 709 


n 535 


n 


2 73 


1 "lO 






794 


n 1 0'^ 


U.'4:UO 


1 8'^p-94 


4 04 






7C 


0.00479 


0.879 


0.0491 


4.15e-20 






/7n 7 ^ kFIP 


3C 




n 61 5 




n 


1 98 


1 "in 


(JOOOtiC/Z/ 


fif! 




n 1 09 






9 Q1 

Zl . X 






7C 


0.00703 


0.412 


0.0689 


6.13e-21 




2 fi 


/In — 90 knr 


3C 


n 0253 


n 802 


n 192 


n 

yj 


2 28 


1 "sn 




fif! 


n 564 


n 021 2 

yj .yj ^ J. £j 


n 533 


1 83p-94 


3 64 






7C 


1.34e-04 


0.893 


0.00708 


8.39e-22 




2 fi 


/In — n-t ^ 
PO — HI 


3C 


n nQ2i 


n n2fin 




n 


n Q87 


1 






n 99Q 


n ni 1 7 




1 ^"^p 94 


1 97 






7C 


6.56e-04 


0.412 


0.0108 


4.30e-20 




9 fi 


f 

/In — rin 






n 495 




n 


1 58 

X .OO 






fiC 


n 2Q3 


7 5np-n4 


n 55fi 


1 8'^p-24 


9 1"^ 

z< . xo 






7C 


7.74e-05 


0.204 


0.00291 


6.13e-21 




9 fi 


irmn[ris) ' 






n nfiQ7 




n 


1 '^1 

X .O X 




OOOOt: C/Zi 


fif! 


n 994 


n 1 09 


n 55fi 




1 9"^ 








z.ooe-Uo 


u. /yu 


n ni no 


o.oye-zz 




2.6 


7mm(/is) = 100 


3C 


8.29e-04 


2.01e-07 


0.00485 





0.199 


150 


3888492 


6C 


0.0809 


4.05e-06 


0.199 


6.76e-22 


0.216 






7C 


9.50e-10 


0.0267 


9.25e-07 


1.74e-18 




2.6 


7maj;(;is) = 3 X 10^ 


3C 


0.0674 


0.212 


0.360 





1.69 


150 


3888492 


6C 


0.869 


0.159 


0.503 


1.83e-24 


2.03 






7C 


0.00307 


0.0465 


0.0491 


6.13e-21 
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Table 2 — Continued 



X 

Tuax^ 


Model 
Ensemble Size 


Survey 


V{P) 


V{D) 


V{z) 


V{a) 






^ — Zi.UU± 




0fi9Q 


001 59 


483 


9 4'ip-39 


1 99 

X . ZZ 








n 

u.u / u 


9 ^^p-Oi^ 
z .tj(j*3 uu 


u.ouz 




1 31 






7C 


1.36e-04 


0.131 


0.00458 


6.90e-17 






n — 9 3 

jJ — z. o 




03^6 


510 


1 35 


n 


9 09 

Z . UZ 


1 1^0 


(JOOOtiC/Z 


fif! 


n 9Q3 


1 09 

U . X uz 


71 fi 


1 83p-94 

± . OOC Zt: 


3 90 






7C 


0.00108 


0.999 


0.0238 


3.72e-20 






f — C] fiyc; 

c — u.u / u 




n 091 1 


001 59 


997 

u.zz / 


n 


1 07 


1 "lO 


OOOOrtc/Z 


fif! 


564 


001 58 
u.uuxuo 


408 


1 83p-94 

± .OOC Zrt 


1 19 






7C 


6.64e-04 


0.190 


0.0161 


6.13e-21 




9 fi 


p — 1 4 




1 99 

\j . J. zz 


filQ 


5fi8 


n 

U 


9 99 

Z . ZZ 


1 


'^RRR4Q9 

OOOO^t/Z 




n 47Q 






1 "^^p 9"^ 

1 .OOt: ZO 


18 






7C 


2.44e-04 


0.199 


0.0108 


4.30e-20 




9 

z, .u 


n — n 9 
/y — u.z 




n nnQ99 


n (T^fi^ 




n 

u 


1 ni 




OOOOtiC/Z 


fif! 


333 


7 50p-04 


533 


1 83p-94 

± .OOC Zrt 


1 17 






7C 


0.0173 


0.204 


0.293 


4.30e-20 




9 


T) — n R 




u.u -L 1 




n 1 Q9 


n 


1 44 


1 "in 

J-UU 


OOOOtiC/Z 


fin 


41 fi 


001 59 


55fi 


1 83p-94 


9 05 
z.uu 








i.iOe-0/ 


n 1 oo 


D.oie-04 


/( 1 ^ on 

4.ioe-zU 




2.6 


r = 2 X 10"^ 


3C 


0.0253 


4.88e-05 


0.112 





0.680 


150 


3888492 


6C 


0.564 


1.65e-07 


0.267 


1.40e-23 


0.696 






7C 


3.58e-06 


0.0273 


4.06e-04 


1.74e-18 




2.6 


r = 3 X 10-=^ 


3C 


0.0497 


0.167 


0.135 





1.37 


150 


3888492 


6C 


0.0959 


0.240 


0.564 


1.83e-24 


2.14 






7C 


1.34e-04 


0.732 


0.00708 


5.91e-21 
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Table 2 — Continued 



X Model 




'P[P,D,z,a] 


TMax^ Ensemble Size Survey V{P) V{D) V{z) 


Via) 


'P[P,2D,z,a] 



2.6 


e = 1.4 


3C 


0.0906 


0.178 


0.580 





1.98 


150 


/3 = 1.6 


6C 


0.416 


0.242 


0.408 


1.40e-23 


2.76 




3888492 


7C 


0.00110 


0.732 


0.0108 


2.79e-19 




2.6 


e = 1.4 


3C 


0.0910 


0.216 


0.581 





2.04 


150 


ao = 7.5 kpc 


6C 


0.340 


0.485 


0.408 


1.40e-23 


2.90 




3888492 


7C 


0.00108 


0.565 


0.0238 


1.74e-18 




2.6 


(3 = 1.6 


3C 


0.207 


0.216 


0.535 





2.26 


150 


flo = 7.5 kpc 


6C 


0.724 


0.340 


0.194 


1.35e-23 


3.13 




3888492 


7C 


0.0106 


0.732 


0.0491 


4.30e-20 





^TMax in units of Myr. 



^Al l other parameters are as in the MMK mod el (^: following iMa nolako u fc Kirk 
(I2OO2I ). with hotspot size increasing according to lJeyakumar fc Saikial (i2000l ). 

^One cannot be confident of the validity of the K-S statistic as the detected sample 
in the simulation is smaller (or, much smaller) than in the actual catalog. 

■^Parameters defining the external environment density profile are set to those of the 
KDA model: /3 = 1.9, Oq = 2 kpc, po = 7.2 x 10-^2 kg m'^. 

Vl = Po (Default)/1.5 = 1.133 X lO^^S ^-3_ 

V2 = 2 X Po (Default) = 3.4 X lO^^s kg m-^. 
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Table 3. K2000 Model: 1-D K-S Statistics for a Few Parameter Variations 



X 


Model 












'P[P,D,z,a] 


a 

Max 


Ensemble Size 


Survey 


V{P) 


ViD) 


V{z) 


Via) 


'P[P,2D,z,a] 


2.6 


Default ^ 


3C 


5.95e-ll 


4.06e-12 


8.21e-07 


3.89e-42 


0.0127 


500 


114900 


6C 


0.00110 


3.59e-09 


0.0185 


1.80e-24 


0.0130 






7C 


9.89e-17 


5.39e-04 


1.24e-09 


1.91e-13 




3.0 


Default 


3C 


8.30e-09 


0.0151 


2.69e-07 


3.89e-42 


0.152 


150 


4861474 


6C 


3.05e-06 


2.13e-04 


0.00264 


1.80e-24 


0.302 






7C 


1.04e-19 


0.221 


1.48e-13 


1.91e-13 




3.0 


Default 


3C 


8.30e-09 


0.0200 


3.53e-06 


2.00e-41 


0.241 


150 


111072 


6C 


8.19e-05 


1.52e-04 


0.00211 


4.92e-21 


0.481 






7C 


4.28e-21 


0.360 


9.20e-14 


2.78e-13 





^Tuax in units of Myr. 

'^AU other dynamical and lobe pow er evolution model parameters are same as in the 
fiducial model of K2000 jKaiseilboOol ). 



Table 4. Model 1-D K-S Statistics using 



iGrimes. Rawlings. fc Willott 



(120041 ) RLF 



X 

rp b 

J- Max 


ModeP 
Ensemble Size 


Survey 


ViP) 


ViD) 


Viz) 


Via) 


V[p^D,z,a] 
V[p^2D,z,a] 


2.6 


KDA 


3C 


5.59e-08 


8.84e-06 


l.Ole-04 


3.01e-08 


0.494 


500 


3712083 


6C 


0.140 


0.00193 


0.310 


3.66e-10 


0.690 










n Q 1 o 
U.oiy 


n nnono 

U.UUZUZ 


n ni nfi 
U.UiUD 




3.0 


KDA 


3C 


0.0281 


0.0976 


0.0583 


3.34e-09 


0.410 


100 


1958652 


6C 


0.107 


0.0204 


0.0397 


1.96e-06 


0.633 








n nnQnQ 
U.UUdUd 


n 1 on 


U.UUoyy 


n nnoQn 




3.0 


KDA 


3C 


0.0611 


0.217 


0.0836 


1.28e-ll 


1.14 


150 


2747159 


6C 


0.394 


0.349 


0.0735 


6.89e-07 


1.82 








U.Ui 1 o 


U./D4 


n ni 1 /I 
U.Uii4 


Q nc^ ne; 

o./oe-Uo 




3.0 


KDA 


3C 


0.133 


0.0125 


0.180 


1.36e-10 


1.44 


200 


3212793 


6C 


0.416 


0.00615 


0.191 


6.89e-07 


2.14 








U.Udoz 


n VQO 
U. 1 oZ 


n nnvofi 
U.UU 1 /o 


n ni nfi 
U.Uiao 




3.0 


KDA 


3C 


0.175 


0.0253 


0.387 


5.04e-14 


1.39 


300 


3813260 


6C 


0.115 


0.349 


0.0245 


6.89e-07 


2.12 








U.UU ( zo 


U.0Z4 


n nn/1 fi/i 
U.UU404 


n nnc^nn 
U.UUoUU 




3.0 


KDA 


3C 


0.122 


0.00101 


0.192 


1.41e-09 


1.32 


500 


5105485 


6C 


0.0454 


0.0968 


0.0255 


1.96e-06 


2.26 






7C 


0.00466 


0.942 


0.0165 


0.00468 




2.6 


BRW 


3C 


4.63e-07 


5.54e-08 


3.43e-05 


8.34e-04 


0.467 


500 


3712083 


6C 


0.0770 


0.0265 


0.212 


3.66e-10 


0.751 






7C 


5.04e-05 


0.438 


9.03e-05 


0.00304 




3.0 


BRW 


3C 


5.84e-06 


2.27e-21 


2.04e-05 


0.00483 


0.303 


100 


1958652 


6C 


0.177 


0.00158 


0.194 


6.76e-07 


0.363 






7C 


6.84e-06 


0.0643 


1.27e-05 


0.0104 
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Table 4 — Continued 



X 


ModeP 
Ensemble Size 


Survey 


V{P) 


V{D) 


V{z) 


V{a) 


V[P,D,z,a] 
V{P,2D,z,a] 




RRW 

UIVV V 


3r! 


3 fi3p-05 
o.uoc uu 


5 73p-1 9 


6 18p-05 

u. XOC^UO 


0398 


fi99 




97471 "iQ 




n 1 9Q 

U. XZiC/ 


u.uuxuo 


9R7 




Q^Q 






7C 


6.67e-04 


0.353 


1.36e-04 


0.0104 






RRW 


3r! 

Ow 


O .Uc/C UO 


c; DOp-I 8 

O .UUL/- X O 




n 0984 

U. UZUt: 


R9R 
u .ozo 




391 97Q3 

OZ/±Z/ / C/O 


U 


n 1 77 

U . ± 1 1 


U .UUUZjZj 


n 987 

U. Z.O ( 


1 8qp-0fi 


1 33 

X .oo 






7C 


2.29e-04 


0.533 


3.99e-04 


0.0194 




o.u 


RRW 

UIW V 


3C 

OV7 


9 Q1p-04 


3 53p-1fi 
o . u oc ± u 


1 05p-04 


0857 

u.uou / 


769 

U. / \JiJ 


ouu 


3S1 39fi0 

OOiOZiUU 




113 

u. ± ±o 


7 50p-04 

/ .out? Ut: 


987 


fi 76p-07 


1 17 

X . X / 






7C 


0.0159 


0.434 


0.00180 


0.0198 




o.u 


RRW 

J_)1VV V 


3r! 

Ov^ 


3 nnp-n4 

O .UUt- Ut: 


9 nnp-91 

Z( . UU" Z X 


9 Q4p-n4 

. C/TiC Ut: 


4 93p-04 

t:. ZOC Ut: 


fiq4 

U .Uc/t: 




u ±uu^ou 


Uv^ 


U .Z,c/0 


(SOp 0^ 

U .UUL.-UO 


9R7 


9 "^"^p 07 


1 09 

X .UZi 






7C 


7.74e-05 


0.353 


7.67e-05 


0.00500 




9 

z, .u 


MK 




^ 89p 1 1 


1 97p DQ 


1 q^p 


n 

u 


719 

U. 1 XZi 


uuu 


371 9083 

O / XZiUOO 




1 fil 
u. xu± 


00903 

U.UUZiUO 


439 


1 80p-94 


1 O'n 
X .uu 






7C 


1.55e-05 


0.548 


0.00126 


1.58e-15 




o.u 


MK 




n ni 99 


n m 87 


n 148 


n 


U .U 1 ^ 


1 on 

xuu 










n 0149 


1 4np-93 


R1 Q 

U.OX t? 








U.UU/4/ 


O.OoZO 


n nm on 
O.OOioy 


z.oye-io 




3.0 


MK 


3C 


0.0494 


0.125 


0.458 





1.36 


150 


2747159 


6C 


0.396 


0.336 


0.0142 


1.40e-23 


1.95 






7C 


0.00181 


0.278 


0.00484 


2.16e-15 




3.0 


MK 


3C 


0.00557 


0.0253 


0.222 





1.26 


200 


3212793 


6C 


0.798 


0.323 


0.0735 


1.02e-22 


1.73 






7C 


0.00185 


0.272 


0.00471 


2.09e-15 
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Table 4 — Continued 



X 


ModeP 
Ensemble Size 


Survey 


V{P) 


V{D) 


V{z) 


V{a) 


V[P,D,z,a] 
V{P,2D,z,a] 




MK 

iVJ-Xv 


3r! 


0950 


00947 


947 


n 

yj 


1 1 1 


ouu 










u.uxou 




1 fil 






7C 


0.0255 


0.422 


0.00479 


2.03e-15 






MK 








474 


n 


8q7 


■inn 

0\j\j 


(J ±VJtjT:0<J 


u 


f) 0441 




1 fiQp-04 


q 8qp-9'? 


1 94 






7C 


0.0484 


0.327 


0.00471 


3.94e-16 






MRRW 


3C 


Q 04p-05 


8 84p-06 


00988 


1 59p-08 


50594 


uuu 


371 9083 




5fi4 


9 55p-05 


18Q 


7 69p-08 


50597 






7C 


0.0170 


5.25e-06 


0.00703 


0.0194 




(J.VJ 


MRRW 




044 


018fi 


0'?90 




4'^5 


1 on 






n 97R 






7 R9p or 








7C 


0.00486 


0.0300 


6.70e-04 


1.41e-05 




o .u 


MRRW 

iVlJ_) JTLV V 








u.uuxu 


1 q^p IzL 


n (SI Q 

U .U X c/ 




97471 5Q 




118 


1 09 


0955 


7 69p-08 








7C 


0.00307 


0.264 


0.00291 


3.32e-05 






MRRW 

1VJ.JJ J. IjV V 




n 1 Q4 


4 87p-0fi 


1Q9 


1 98p-1'^ 

X . Z( O " X (J 


q'^8 








n 797 


01 0'^ 




7 fi9p-08 


950 








0.00450 


O.OOoyo 


0.010/ 


9.01e-0o 




3.0 


MBRW 


3C 


0.253 


0.00159 


0.192 


5.09e-12 


0.707 


300 


3813260 


6C 


0.276 


0.0117 


0.0735 


7.62e-08 


0.734 






7C 


0.00481 


0.0194 


0.0111 


0.00500 




3.0 


MBRW 


3C 


0.385 


1.42e-06 


0.312 


9.59e-07 


0.882 


500 


5105485 


6C 


0.120 


0.00150 


0.0255 


7.62e-08 


0.923 






7C 


0.0358 


0.0430 


0.0160 


0.00976 
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Table 4 — Continued 



X 


Model'^ 












'P[P,D,z,a] 


rp b 

Max 


Ensemble Size 


Survey 


V{P) 


V{D) 


V{z) 


V{a) 


'P[P,2D,z,a] 


2.6 


MMK 


3C 


0.0681 


0.0125 


0.474 





1.15 


500 


3712083 


6C 


0.724 


0.00300 


0.177 


1.83e-24 


1.20 






7C 


0.00451 


0.0294 


0.00710 


4.15e-20 




3.0 


MMK 


3C 


0.0200 


0.00851 


0.00903 





0.210 


300 


3813260 


6C 


0.107 


0.154 


0.00189 


1.83e-24 


0.316 






7C 


0.00750 


0.00102 


0.00114 


4.00e-20 




3.0 


MMK 


3C 


0.103 


0.0696 


0.0206 





0.219 


500 


5105485 


6C 


0.00775 


1.44e-04 


1.19e-05 


1.83e-24 


0.298 






7C 


0.0115 


0.00920 


0.00188 


5.91e-21 





^AU the dynamical and radio lobe power evolution model parameters arc same as in the 
default version of the corresponding model (KDA, BRW, MK, MBRW, MMK). 

^TMax in units of Myr. 
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Table 5. 2-D K-S Test Results for the Two Modified Models 



Model 
Parameters 


Survey 




2-D KS Probability, P(K-S) 

V{P-D) V{z-D) V{P-a) V{z-a) 


V{D-a) 


MBRW 


3C 


0.00234 


4.17e-10 


2.10e-08 


8.27e-08 


6.07e-08 


5.47e-12 


Default 


6C 


0.458 


0.00506 


0.00667 


1.46e-06 


1.48e-06 


1.53e-04 




7C 


0.00309 


8.54e-04 


0.00222 


0.00733 


0.00648 


0.0112 


MBRW 


3C 


0.299 


0.00119 


0.00246 


5.99e-ll 


3.86e-10 


4.64e-14 


Varied ^ 


6C 


0.781 


0.434 


0.496 


3.21e-06 


7.85e-06 


2.76e-04 




7C 


1.60e-04 


0.00297 


0.00237 


6.00e-05 


6.28e-05 


0.00228 


MMK 


3C 


0.0366 


1.48e-04 


0.00444 


2.29e-40 


2.28e-39 


1.88e-32 


Default ^ 


6C 


0.0656 


0.0325 


0.183 


1.70e-17 


2.32e-18 


3.76e-16 




7C 


0.0151 


0.00180 


0.00355 


3.53e-15 


5.25e-14 


2.88e-13 


MMK 


3C 


0.147 


0.143 


0.421 


8.35e-40 


5.53e-37 


8.64e-30 


Varied 


6C 


0.254 


0.282 


0.235 


1.53e-17 


2.08e-18 


2.22e-14 




7C 


0.0203 


0.00627 


0.0290 


2.27e-14 


2.96e-13 


2.91e-12 



^^Simulations with the respective parameters of MBRW model as used in §21 Initial 
population generated using x = 2.6, Tmux = 500 Myr. 

"^MBRW model simulation using initial population with x = 3.0, TMax = 300 Myr. The 
power evolution is with parameter change tbf = 100 yr, other parameters set to their default 
values, for a case with initial source population size = 4963343 (the last but 6th entry in 
Table [I]). 

'^Simulations with the respective parameters of MMK model as used in ^ Initial popu- 
lation generated using x = 2.6, Tmux = 500 Myr. 

■^MMK model simulation using initial population with x = 2.6, Tjiiax = 150 Myr. The 
power evolution is with parameter change (3 = 1.6, other parameters set to their default 
values, for a case with initial source population size = 3888492 (9th entry in Table [2]). 
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Table 6. 4- variable Spearman Partial Rank Correlation Analysis 



Data Model (combining all surveys ) 

MBRW MMK 



Coeff. 


All " 


Default 


Varied ^ 


Default 


Varied ^ 


T" D n 


0303 


0.183 


0309 


0.160 


0.127 


^PD,za 


0.478 


2.94 


0.489 


2.57 


2.02 




0.716 


0.667 


0.649 


0.0754 


-0.103 


^Pz,Da 


14.2 


12.8 


12.2 


1.20 


-1.64 


fDz,Pa 


-0.268 


-0.348 


-0.237 


0.218 


0.00189 


'^Dz,Pa 


-4.33 


-5.77 


-3.83 


3.52 


0.0300 


f Pa,Dz 


0.147 


-0.0163 


-0.0139 


-0.640 


-0.750 


'^Pa,Dz 


2.33 


-0.259 


-0.219 


-12.1 


-15.4 


fDa,Pz 


0.472 


0.0649 


-0.275 


0.605 


0.346 


^Da,Pz 


8.08 


1.03 


-4.47 


11.2 


5.72 


fza,PD 


-0.0234 


0.132 


0.194 


-0.590 


-0.609 


'^za,PD 


-0.369 


2.12 


3.12 


-10.8 


-11.2 



^The four observables P, D, z and a for the 3C, 6C and 7C 
III surveys (whether real or simulated), combined together in 
a single sample. 

^The particular parameters used are the same as those in 
Table \5\ for each of the MBRW and MMK varied models. 

'^Spearman partial rank correlation coefficient between two 
variables P and D, when the other two variables z and a are 
kept fixed. 
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Significance level associated with the correlation between 
P and D, independent of z and a. 
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Table 7. Relevant Volume Fractions for Selected Models and Modifications 



X 




Ensemble 


Model 


Ratioac A (xlO'^) 


Volume 


Total 




(Myr) 


Size 


Parameters 




(Mpc-3) 


Fraction, l 


Vol-Frac, ( 










BRW model 








2.6 


250 


1466378 


Default 


2.12 


0.931 


0.00147 


0.00623 


2.6 


500 


1561417 


Default 


1.15 


0.532 


0.0123 


0.0301 





950 




Default 


1.33 


0.997 


n no904 




O.yJ 






flo = 7.5 kpc 


0.869 


2.13 


n nnQfi4 

U. UUC/Ut: 


040S 


O.yJ 


oyjyj 




Default 


1.21 


0.997 


n m 7Q 










Oo = 7.5 kpc 


1.03 


2.19 


0667 

w . www 1 


W . X w»_> 








MBRW model 








2.6 


500 


4397469 


Default 


2.87 


1.50 


0.0139 


0.0339 


3.0 


300 


4963343 


hf = 100 yr 


b 1.94 


2.66 


0.00658 


0.0240 










KDA model 














Default 


2.17 


1.60 


n noo'^94 


0091 5 






4SQ74fiQ 


Default 


2.18 


1.50 


01 83 

W . W X (JO 


0459 


3 


150 

X (J w 


1618248 


Default 


0.703 


1.66 


000796 

W . W W W 1 t^W 


00530 

W . W Wt-.^OW 


o.u 


1 ^0 




P, Po "''^ 


0.993 


A QQ 


u. uuoo ± 


n 0990 


3.0 


500 


3419466 


Default 


0.559 


1.16 


0.0428 


0.106 


3.0 


500 


6451283 




0.393 


2.19 


0.225 


0.555 










MK model 








2.6 


150 


3888492 


Default 


3.26 


3.99 


0.000772 


0.00512 


2.6 


500 


4397469 


Default 


1.37 


1.50 


0.0286 


0.0699 


3.0 


150 


4861474 


Default 


1.20 


4.99 


0.00213 


0.0141 


3.0 


150 


4861474 


'~fmax{hs) 


1.11 


4.99 


0.00230 


0.0153 


3.0 


500 


4886474 


Default 


0.421 


1.66 


0.0844 


0.206 


3.0 


500 


4886474 


lfmax{hs) 


0.414 


1.66 


0.0858 


0.209 










MMK model 








2.6 


150 


3888492 


/3= 1.6 s 


1.59 


3.99 


0.00220 


0.0146 
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Table 7— Continued 



^ ^Max 


Ensemble 


Model 


Ratio3c 


A, (xlO-S) 


Volume 


Total 


(Myr) 


Size 


Parameters 




(Mpc-3) 


Fraction, l 


Vol-Frac, ( 


2.6 500 


4397469 


Default 


0.710 


1.50 


0.0552 


0.135 



^1-D K-S best-fit case of BRW. 

K-S best-fit case of MBRW. 
'^Parameter variations: p — 2.12 and po = Po (Defauit)/2 = 3.6 x 10"^^ kg m 

K-S best-fit case of KDA. 
^'Parameter variation: 'jrnax(hs) = 3 x 10^. 
^l-D K-S best-fit case of MK. 
si-D K-S best-fit case of MMK. 



